Coherent control of interacting electrons in quantum dots by traveling the spectrum 
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Quantum control of the wave function of two interacting electrons confined in quasi-one- 
dimensional double- well semiconductor structures is demonstrated. The control strategies are based 
on the knowledge of the energy spectrum as a function of an external uniform electric field. When 
two low-lying levels have avoided crossings our system behaves dynamically to a large extent as a 
two-level system. This characteristic is exploited to implement coherent control strategies based on 
slow (adiabatic passage) and rapid (diabatic Landau-Zener transition) changes of the external field. 
We apply this method to reach desired target states that lie far in the spectrum from the initial 
state. 

PACS numbers: 73.63.-b, 78.67.Hc 
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The control of quantum systems is a fundamental chal- 
lenge in physical chemistry, nanoscience, and quantum 
information processing [l|, Quantum control is the 
manipulation of the temporal evolution of a system in 
order to obtain a desired target state or value of a cer- 
tain physical observable. From the experimental point 
of view, the techniques of quantum control are highly 
developed in the area of magnetic resonance, and more 
recently great progress has been made in quantum chem- 
istry thanks to the development of ultrafast laser pulses 
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Coherent control in semiconductor quantum dots has 
become an active field of research in the last 15 years. 
Early works on electron localization in double well sys- 
tems spurred intense theoretical activity. In a seminal pa- 
per, Grossmann et al. [3| showed that, by applying an ap- 
propriate AC electric field, the tunneling of the electron 
between the wells could be coherently destroyed, thereby 
maintaining an existing localization in one of the wells. 
Shortly after, Bavli and Metiu Q found ways to, starting 
from the delocalized ground state, localize the electron 
wave function and then to preserve the localization with 
a precisely taylored time-dependent electric field. A large 
body of literature followed these pioneering works. A 
decade later, the first steps in the theoretical exploration 
of localization and control of two interacting electrons in 
quantum dots were made [(| 0, 0] • Whereas Zhang and 
Zhao studied a model two-level system, Tamborenea and 
Metiu studied a more realistic multi-level system inspired 
by quasi-one-dimensional semiconductor nanorods. The 
study of two-electron localization and control in double 
dots has remained active ever since @, [13, 11, 12, HI- 

In this Letter we propose an efficient method to control 
the wave function of two interacting electrons confined 
in quasi-one-dimensional nanorods [g, [ijj]. The control 
method is based on the knowledge of the energy spec- 
trum as a function of an external uniform electric field. 
The method requires that the system behaves locally — 
near avoided level crossings — as the Landau-Zener (LZ) 



two- level model [l^]. This fact is exploited to navigate 
the spectrum using slow (adiabatic) and rapid (diabatic) 
changes of the external control parameter. Although this 
characteristic may seem rather restrictive, it is, in fact, a 
general property of systems with interaction between its 
energy levels. The level repulsion must not be too strong, 
though, so that the spacing at the avoided crossings re- 
mains smaller than the mean level spacing. 

Let us consider a quasi-one-dimensional double quan- 
tum dot with two interacting electrons in it in the pres- 
ence of a spatially uniform electric field [8( . Denoting by 
z the longitudinal coordinate, the Hamiltonian of the two 
electrons reads 



H = — 



2m \dz 

VC{\Z! - 



*2|) 



-V{z x ) + V{z 2 ) 
e(z! + z 2 )E(t), 



(1) 



where m is the effective electron mass in the semicon- 
ductor material, Vc is the Coulomb interaction between 
the electrons, V is the confining potential, and E(t) is 
an external time-dependent electric field. The confining 
potential is a double quantum well with well width of 
28 nm, interwell barrier of 4 nm, and 220 meV deep (a 
typical depth for a GaAs-AlGaAs quantum well). We 
mention that the control results that we report here are 
robust with respect to the fine tuning of the parameters 
of the structure. We normally assume that the initial 
state is the ground state, which is a singlet. Since the 
Hamiltonian is spin independent, the total spin is con- 
served and the spatial wave function remains symmetric 
under particle exchange at all times. 

We first consider the case of constant electric field, to 
be considered at this point as a parameter in the Hamil- 
tonian. In Fig. [1] we show the spectrum of eigenener- 
gies of the Hamiltonian of Eq. (Q]) versus electric field. 
The energies ti(E) and eigenstates 4>i(E, z x , z 2 ) are ob- 
tained by numerical diagonalization. We have used as 
basis set the symmetric combinations of the twelve low- 
est single-particle eigenf unctions. Thus, our basis of the 
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two-particle Hilbert space has 12*(12+l)/2=78 states. 
We can see in Fig. Q] that the spectrum is composed 
of nearly straight lines 16]. A closer look reveals that 
all level crossings are avoided, giving rise to adiabatic 
curves that never cross each other. The fact that all 
crossings are avoided ones is due to the electron-electron 
interaction, which couples the energy levels of the non- 
interacting system [l7j |. 

The straight lines of the spectrum are distributed in 
three clearly distinguishable groups: those with negative, 
zero and positive slope. In each group the slopes are 
very similar, and far from the avoided crossings, the wave 
functions have a distinct kind of localization, as follows 
(see Figs. Q] and HJa-f)): 

i) for zero slope the electrons are delocalized, that is, each 
electron is in a different dot (Fig. [2](a,b)), 

ii) for negative slopes both electrons are in the left dot 
(Fig.dcd)), and 

iii) for positive slopes both electrons are in the right dot 
(Fig. EM)). 

Along a given straight line the eigenstates do not change 
much, thus each line has associated a characteristic shape 
of the wave function. Near an avoided crossing, states 
with different types of localization mix. These mixed 
states do not have a well-defined localization type. As 
an example, in Fig. g) and (h) we show the eigenstates 
at the avoided crossing between the first two levels at 
E = 4.81 kV/cm and between the levels 21 and 22 at 
E = 8.48 kV/cm, respectively. 
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FIG. 1: Energy spectrum of our system as a function of the 
external uniform electric field. 

Our goal is to find a method to control the wave func- 
tion. The previous discussion about the spectrum and 
the characteristics of the wave functions at and far from 
the avoided crossings suggests a possible control strategy. 
For example, starting at state (a) of Fig. [T] and varying 
slowly (adiabatically) the electric field we reach state (c), 



which has a different type of localization (see Fig. [2]). On 
the other hand, if we vary the parameter E quickly the 
final state will have the same localization as the initial 
one (see below, Fig. [3]). These types of transitions will 
be the building blocks of our control strategy. The exact 
meaning of slow and fast in this context is given by the 
LZ model 0. 

The LZ theory treats a rather simple and generic situ- 
ation of two-level avoided crossing. The LZ Hamiltonian 
in the diabatic basis {|+), |— }} is 



H = 



e+(A) e 
e e_(A) 



(2) 



where e is a constant while the diagonal £±(A) = Eq + 
a±\ are linear functions of the parameter A. It is as- 
sumed that the parameter A varies linearly with time. It 
can be shown that the probability of remaining in state 
|+) is given by 



exp 



-2TT£ 2 



a-)\ 



(3) 



where A is the rate of variation of the parameter. Con- 
versely, P- = 1 — P+. By selecting the rate of change 
of the parameter A we can control the final state of the 
system. For slow variations of A (A <C %-)) the system 
follows the adiabatic curve passing from the initial dia- 
batic state to the other one. In the opposite case, when 

2 

A> j, the evolution takes place on the diabatic curve 
and the state remains as the initial one. It can be shown 
that near an avoided cros sing our system can be treated 
as a two-level LZ model [18(. We compared the evolu- 
tion of the two- level LZ Hamiltonian (Eq. |[2}) to that 
of the complete Hamiltonian (Eq. |T])) for many rates of 
variation of the parameter (electric field) and we found 
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FIG. 2: Spatial wave functions <f>i{E, zi, z<i) corresponding to 
labels (a) to (h) of Fig. [T] States (a) and (b) have one electron 
in each well, wave functions (c) and (d) are localized in the 
left well, and wave functions (e) and (f ) have both electrons in 
the right well. States (a) to (f) are far from avoided crossings 
and therefore have well-defined localization properties. This 
is not the case for eigenstates (g) and (h) , which are at avoided 
crossings. 
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very close agreement. The numerical solution of the time- 
dependent Schrodinger equation was obtained using the 
usual fourth-order Runge-Kutta method [f| Q • We use 
a time step of 0.1 fs, which ensures that the precision of 
all the reported probabilities and overlaps is better than 
0.4%. 

The adiabatic transition from state (a) to (c) (see Figs. 
Q] and [2] (a,c)) is an extremely simple solution to the 
problem of localization in a realistic double quantum dot 
@j B 0- In fact, more complex and interesting control 
problems can be tackled. Namely, by suitable combina- 
tions of rapid (diabatic) and slow (adiabatic) changes of 
the electric field it is possible to travel over the spec- 
trum connecting distant pairs of states. That is, we can 
not only control the distribution of electrons in the dots 
(localization and delocalization), but we are also able to 
choose finer details of the final state, like, for example, 
its nodal domains. 

We now show a complex navigation through the spec- 
trum in which the system is taken from the ground state 
to a specific high-energy state. In this case, the target is 
state (b) of Fig. [2] (see also Fig. [I]), which shares the local- 
ization type with the ground state (they are both delocal- 
ized), but has a more complex nodal structure. This pro- 
cess is depicted in Fig. [3] We want to reach the excited 
state by means of diabatic and adiabatic transitions. The 
intended navigation path is displayed with arrows in the 
spectrum of Fig. [3^a) . The small arrows indicate slow 
variations of the control electric field whose objective is 
to follow the adiabatic states (i.e. the eigenfunctions of 
the Hamiltonian at the successive values of the electric 
field) . The long arrows denote diabatic transitions at the 
avoided crossings. Here we use instantaneous jumps of 
the electric field, but we have checked that risetimes of 
the order of 0.1 ps do not change significantly our results. 
The time dependence of the electric field is shown in Fig. 
[3](b). In order to know in detail how the evolution of 
the wave function proceeds, we show in Fig.[3jc,d) differ- 
ent aspects of the time-dependent wave function. In Fig. 
EHc) we compute the time-dependent localization proba- 
bilities, Pll, that the two electrons are in the left well, 
and Prl, that they are in different wells. Fig. [3jd) gives 
the absolute value of the overlap of the evolving state 
with the energy eigenstate <j>i(E), and Fig. [3fe) displays 
the spatial wavefunction at the various times indicated 
in (d). 

Let us discuss in detail the behavior of the evolving 
wave function. We start from the ground state with zero 
electric field and move adiabatically up to a field of 3.5 
kV/cm, as shown in Fig. [3]Ja). This process takes 100 
ps (see Fig. EJb)). Note in Fig. [3jd) that the overlap 
with the adiabatic eigenstate </>i is approximately 1. Af- 
ter that the field is quickly increased (the first jump of 
the electric field in Fig. [U(b), at t = 100 ps) so that the 
system goes through the first avoided crossing. At this 
avoided crossing the overlap between the evolving state 



and the adiabatic eigenstates begins to decay, from al- 
most 1 to close to 0.97. For example, at time t\ = 190 
ps the wave function still resembles the ground state (see 
Fig. EJe)) but displays a degree of mixture with localized 
states, as evidenced by the finite probability in the lower- 
left quadrant. Afterwards, we move slowly until = 300 
ps passing adiabatically the second avoided crossing at 
E = 11.2 kV/cm. The mixed nature of the wave func- 
tion at the avoided crossing (at t 2 — 243 ps) can be seen 
in Fig. [U(e). After the crossing, at ts, the wave function 
is highly localized on the left dot (see Fig.[3fe)). At this 
point, we still have ahead of us a long way to the desired 
final state, most of it along a spectral "line" of negative 
slope. Although it would be tempting to make a sudden 
change of the control parameter to its final value, that 
strategy is not satisfactory. It is best to proceed slowly 
far from the avoided crossings and rapidly around them. 
This procedure has the advantage of allowing our state 
to adjust to the gradual changes of the adiabatic states in 
the regions between crossings. These gradual changes do 
not entail changes in the type of localization, but rather 
displacements of the probability density within the quad- 
rants which are already populated. In the remainder of 
the path we cross three avoided crossings diabatically and 
the last one slowly. The latter, being very narrow, turned 
out to be the hardest one to pass adiabatically. This pas- 
sage took 300 ps, as seen in Fig. \3[b). In this avoided 
crossing the wave function goes from a strong localiza- 
tion in the left dot (wave function at £4 of Fig. 02e)) to 
a clear delocalization in the final state. We remark that 
the overlap of the final state with the desired state, at 
< 6 = 800 ps, is 0.93 (see Fig.^d)) 

Before concluding, it is worth comparing the present 
control method with others found in the literature. In 
Ref. [H, a different localization scheme was proposed for 
the system studied here. In that work the electron lo- 
calization is obtained, starting from the ground state at 
zero electric field, with a piecewise constant electric field. 
The time scales involved, of a few picoseconds, which are 
experimentally accesible, are the same in the two ap- 
proaches. However, the present method has the advan- 
tage of being more robust since it does not require the 
fine tuning of the field jumps and, furthermore, it gives 
higher probabilities of reaching the target state. More- 
over, in Ref. [1] the controlled state is a superposition 
of only the two lowest-energy eigenstates, while in our 
proposed method a much more general navigation of the 
spectrum is naturally possible, allowing to connect very 
distant states. More recently, an efficient method of con- 
trol was proposed and applied to the creation of entan- 
gled states in spin systems [l9]. The method is based 
on the navigation of the energy spectrum, expressed as a 
function of an external control parameter, which is var- 
ied always adiabatically. A key element of the controlled 
system of that study is that the interaction, which con- 
trols the nature of the level crossings, can be switched 
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FIG. 3: (a) Schematic diagram of the intended path followed by the state. Short (long) arrows indicate slow (fast) variations of 
the electric field. The initial state is marked as ■ and the target state as ©. (b) External electric field (the control parameter) 
as a function of time, (c) Probabilities that the two electrons are in the same (left) well {Pll) and in different wells (Prl)- 
(d) Absolute value of the overlap of the evolving state with the energy eigenstate <j>i(E). (e) Spatial wavefunctions at various 
times indicated in (d). 



off at will. This additional control tool is clearly not 
generic. Here we propose a control method based on 
similar general ideas but which is more general in the 
sense that the interaction among levels does not need to 
be controlled. Finally, we mention that nonadiabatic LZ 
transitions have been applied to control single-electron 
transport [2(J, and adiabatic passage has been used to 
control localization and suppression of tunneling 21 1 . 

In summary, we have proposed an efficient method to 
control the wave function of a two-electron quantum dot 
system. This method is based in slow and fast varia- 
tions of a control parameter which in our study is given 
by an external electric field that couples to the elec- 
trons through the dipole Hamiltonian. The success of 
this method relies on the condition that the interaction 
between neighboring levels be well described by the two- 
level LZ model. Although this may seem a severe limi- 
tation to the applicability of the method, this condition 
is quite generally satisfied. For example, it has been re- 
cently shown in the paradigmatic stadium billiard that 
the transitions between neighboring levels are of the LZ 
type when the billiard boundary is deformed (22|. We 



have shown that with this method an effective control of 
the wave function can be achieved, even for rather am- 
bitious control goals. While the technology to produce 
the required ramped electric fields, controllable in the 
picosecond time scale with enough flexibility, may not 
be currently available, we expect that our proposal will 
further motivate their experimental development. 
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